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Abstract 

In some communication networks, such as passive RFID systems, the energy used to transfer 
information between a sender and a recipient can be reused for successive communication tasks. In fact, 
from known results in physics, any system that exchanges information via the transfer of given physical 
resources, such as radio waves, particles and qubits, can conceivably reuse, at least part, of the received 
resources. 

This paper aims at illustrating some of the new challenges that arise in the design of communication 
networks in which the signals exchanged by the nodes carry both information and energy. To this end, a 
baseline two-way communication system is considered in which two nodes communicate in an interactive 
fashion. In the system, a node can either send an "on" symbol (or "1"), which costs one unit of energy, 
or an "off signal (or "0"), which does not require any energy expenditure. Upon reception of a "1" 
signal, the recipient node "harvests", with some probability, the energy contained in the signal and stores 
it for future communication tasks. Inner and outer bounds on the achievable rates are derived. Numerical 
results demonstrate the effectiveness of the proposed strategies and illustrate some key design insights. 

Index Terms — Two-way channel, interactive communication, energy transfer, energy harvesting. 

I. Introduction 

The conventional assumption made in the design of communication systems is that the energy 
used to transfer information between a sender and a recipient cannot be reused for future commu- 
nication tasks. There are, however, notable exceptions. An example is given by communication 
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Figure 1. Two-way noiseless binary communication with energy exchange. The total number of energy units is fixed (to five in 
the figure) and transmission of a "1" symbol transfers energy from the sender node to the recipient. See Fig. [3] for a generalized 
model. 



based on wireless energy transfer, such as passive RFID systems HI or some body area networks 
ED, in which a terminal can transfer both information and energy via the transmitted radio signal, 
and the delivered energy can be used for communication by the recipients. For instance, a passive 
RFID tag modulates information by backscattering the radio energy received from the reader 
(see, e.g., 0]|). Another, less conventional, example is that of a biological system in which 
information is communicated via the transmission of particles (see, e.g., O), which can be later 
reused for other communication tasks. A further potential instance of this type of networks is 
one in which communication takes place via the exchange of quantum systems, such as photons, 
which may measured and then reused H|. 

To summarize, any system that exchanges information via the transfer of given physical 
resources, such as radio waves, particles or qubits, can conceivably reuse, at least part, of 
the received resources for later communication tasks. This conclusion is supported by physical 
considerations and practically demonstrated by the existing systems based on this principle 
ifTI , BU. It is emphasized that the possibility to deliver jointly energy and information promises 
not only to ease the energy requirements of various communication systems, but also, more 
importantly, to enable novel applications, such as the body area networks studied in HI, 0. 
Moreover, an understanding of the interplay between energy and information flows could lead 
to insights on the workings of some communication systems in nature 0). 

A. State of the Art 

While the interaction between energy and information continues to be subject of research in 
the physics community (see, e.g., [0), the topic has been tackled from a communication and 
information theoretic level only in a handful of pioneering works, as reviewed in the following. 
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The references 0, flU, [|9]| have focused on the problem of maximizing the information rate 
of a point-to-point system subject to minimum received energy constraints. Specifically, in Q 
a single point-to-point channel was studied, while (8]|, O investigated a set of parallel point- 
to-point channels. To illustrate the trade-offs between the transfer of energy and information in 
a point-to-point channel consider the noiseless transmission of a 4-PAM signal in the alphabet 
{—2,1,1,2}. If one requires the received energy to be the maximum possible, that is, to be 
equal to 4, the maximum transferable information rate is 1 bit per symbol, since one is forced 
to communicate only with the larger energy symbols {—2, 2}. However, with no receive energy 
constraint, one can clearly convey 2 bits per symbol by choosing all available symbols with 
equal probability. This example also explains the substantial difference between the problems 
studied in Q, HI, and that with maximum receive energy constraints studied in IfTOlR The 
optimization of beamforming strategies under a receive energy constraint was tackled in [fTT]] . 
|fT2l for multiantenna broadcast channels. Considerations on the design of the receiver under 
the constraint that, when harvesting energy from the antenna, the receiver is not able to use the 
same signal for information decoding, can be found in lfT3i 

B. Contributions 

In all of the previous work summarized above, the requirement on the energy harvested from 
the received signal is considered to be an additional constraint imposed to the system design. 
This work is instead motivated by the observation that, in more complex network scenarios, 
as mentioned above, the energy harvested from the received signal may be reused for future 
communication tasks. In this case, the energy and information content of the exchanged signals 
should be engineered so as to best suit the requirements of the communication network. To study 
this aspect, we consider a baseline two-way communication system, as illustrated in Fig. [TJ This 
is incidentally the same topology selected by Shannon to initiate the study of networks from an 
information theoretic perspective lfl5l . In the considered model, the two nodes interact for the 
exchange of information and can harvest the received energy. 

To enable analysis and insights, we assume that the two parties involved have a common 
clock and that, at each time, a node can either send an "on" symbol (or "1"), which costs one 

'With a maximal receive energy constraint of 4, one can still clearly transmit 2 bits per symbol. 
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unit of energy, or an "off signal (or "0"), which does not require any energy expenditure. Upon 
reception of a "1" signal, the recipient node can harvest, possibly with some loss, the energy 
contained in the signal and stores it for future communication tasks. In order to introduce the 
main concepts with the minimum of the notation and technical complications, we first consider 
the case in which the two nodes start with a given number of energy units in their batteries, 
which can neither be lost or replenished from outside, and the binary channel in either direction 
is noiseless. 

To see that even this simple scenario offers relevant research challenges, we observe the fol- 
lowing. If there were no limitation on the number of energy units, the nodes could communicate 
1 bit per channel use in either direction given that the channels are ideal. However, if there is, 
say, one energy unit available in the system, only the node that currently possesses the energy 
unit can transmit a "1", whereas the other node is forced to transmit a "0". Therefore, the 
design of the communication strategy at the nodes should aim not only at transferring the most 
information to the counterpart, but also to facilitate energy transfer to enable communication in 
the reverse direction. We study this problem, described in Sec. HI by deriving inner and outer 
bounds on the achievable rate region as a function of the available energy units in Sec. Unl The 
main results are then extended to a model that accounts for energy replenishments and losses, 
along with noisy channels. The generalized model is presented in Sec. [TV] and the generalized 
results are presented in Sec. |V] 

It is finally observed that the class of problems at hand, in which terminals can harvest energy 
from the received signals is related to the increasing body of work on energy harvesting (see, 
e.g., [14J and references therein). However, in this line of work, the energy is assumed to be 
harvested from the environment in a way that is not affected by the communication process, 
unlike the scenario under study. 

Notation: [m, n] = {m, m + 1, n} for integers m < n. N is the set of integer numbers. We 
use the standard notation in ITT61 for information theoretic quantities such as entropy and mutual 
information. If the distribution is Bern(p) we will also write H(p) for the entropy. Capital letters 
denote random variables and the corresponding lowercase quantities denote specific values of 
the random variables. X 1 for an integer i denotes the vector X 1 = JQ). 
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II. System Model 

We consider the binary and noiseless two-way system illustrated in Fig. 1, in which the total 
number of energy units in the system is equal to a finite integer number U > 1 at all times 
and the channels between the two nodes are noiseless. In Sec. [IVl the model will be extended 
to include stochastic stochastic energy losses and replenishments along with noisy channels. At 
any given time instant k, with k E [l,n], the state of the system (Ui t k, U2,k) £ N 2 is given by 
the current energy allocation between the two nodes. Specifically, a state (E/i,*,, U 2 ,k) indicates 
that at the A;th channel use there are Uj^ energy units at Node j, with j = 1,2. Since we assume 
here that Ui ; k + U 2) k = u f° r eacri channel use k E [l,n] (i.e., no energy losses occur), then, in 
this section and in the next, we will refer to as the state of the system, which always imply 
the equality U^k = U — U\ t k- 

At any channel use k E [l,n], each Node j can transmit either symbol Xj ; k = or symbol 
Xj jk = 1, and transmission of a "1" costs one energy unit, while symbol "0" does not require any 
energy expenditure. Therefore, the available transmission alphabet for Node j, j = 1, 2 during 
the kth channel use is 

X U = {0,1} if U hk = u>\ (la) 
and X = {0} if U j>k = 0, (lb) 

so that Xj t k € X u if Uj t k = u energy units are available at Node j. The channel is noiseless so 
that the received signals at channel use k are given by 

Yi, k = X 2 , k and Y 2:k =X^ k (2) 

for Node 1 and Node 2, respectively. 

Transmission of a "1" transfers one energy unit from the sender node to the recipient node. 
Therefore, the state of Node 1 for k E [l,n] evolves as follows 

Ui } k = (Ui^-i — X^k-i) + X 2j fc_i, (3) 

where we set Ux,i = Wi,i < U as some initial state and U 2 ,k = U — U\^- We observe that the 
current state U% t k is a deterministic function of the number U of total energy units, of the initial 
state U^i and of the previously transmitted signals X^ 1 and X^ 1 . We also note that both 
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nodes are clearly aware of the state of the system at each time since U 1:k + U 2)k = U is satisfied 
for each channel use k. 

Node 1 has message Mi, uniformly distributed in the set [l,2 ni?1 ], to communicate to Node 
2, and similarly for the message M 2 £ [1, 2 nR2 ] to be communicated between Node 2 and Node 
1. Parameters Ri and R 2 are the transmission rates in bits per channel use (c.u.) for Node 
1 and for Node 2, respectively. We use the following definitions for an (n, R\, R 2 , U) code. 
Specifically, the code is defined by: the overall number of energy units U; two sequences of 
encoding functions, namely, for Node 1, we have functions f l fc for k e [l,n], which map the 
message M\ and the past received symbols X\~ l (along with the initial state) into the currently 
transmitted signal X^ k e X Ulk ; similarly, for Node 2, we have functions £ 2ik for k e [l,n], 
which map the message M 2 and the past received symbols X^~ l (along with the initial state) 
into the currently transmitted signal X 2yk e Xu 2 k > an ^ two decoding functions, namely, for Node 
1, we have a function gi, which maps all received signals X^ and the local message M\ into 
an estimate M x of message M 2 ; and similarly, for Node 2, we have a function g 2 , which maps 
all received signals X™ and the local message M 2 into an estimate Mi of message Mi. 

We say that rates (R±, R 2 ) are achievable with U energy units if there exists an (n, Ri, R 2 , U) 
code for all sufficiently large n that guarantees reliable communication. We are interested in 
studying the closure of the set of all the rate pairs (R 1 ,R 2 ) that are achievable with U energy 
units, which we refer to as capacity region C(u). Given the noiseless nature of the channels, we 
note that the initial state U^i = u ltl < U does not affect the rate region since in a finite number 
of steps it is always possible to redistribute the energy according to any desired state. 

III. Inner and Outer Bounds 
In this section, we derive inner and outer bounds to the capacity region. 

A. Inner Bounds 

In order to gain insights into the nature of the problem under study, we consider here various 
communication strategies. We start by the simplest, but intuitively important, case with U= 1, 
and we then generalize to U> I. 
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1 ) U= 1 Energy Unit: We start with the special case of one energy unit (U = 1) and assume 
the initial state Uii = 1, so that the energy unit is initially available at Node 1. The other case, 
namely ui \ = 0, can be treated in a symmetric way. In this setting, during each channel use, 
"information" can be transferred only from the node where the energy unit resides towards the 
other node, and not vice versa, since the other node is forced to transmits the "0" symbol. This 
suggests that, when U = 1, the channel is necessarily used in a time-sharing manner, and thus 
the sum-rate is at most one bit per channel use. The first question is whether the sum-rate of 1 
bit/c.u. is achievable, and, if so, which strategy accomplishes this task. 

A Naive Strategy: We start with a rather naive encoding strategy that turns out to be insufficient 
to achieve the upper bound of 1 bit/c.u.. The nodes agree on a frame size F = 2 b > 1 channel 
uses for some integer b and partition the n channel uses in n/F frames (assumed to be an integer 
for simplicity). The node that has the energy unit at the beginning of the frame communicates 
b = log 2 F bits to the other node by placing the energy unit in one specific channel use among 
the F = 2 b of the frame. This process also transfers the energy unit to the other node, and the 
procedure is repeated. The sum-rate achieved by this scheme is 

Ri + R 2 = [bits/c.u.], (4) 

r 

which is rather inefficient: the maximum is achieved with F = 2, leading to a sum-rate of 
+ # 2 = 1/2 bits/ c.u.. 

The previous strategy can be easily improved by noting that the frame can be interrupted after 
the channel use in which the energy unit is used, since the receiving node can still decode the 
transmitted b bits. This strategy corresponds to using a variable-length channel code. Specifically, 
we can assign, without loss of optimality within this class of strategies, the codeword "01" to 
information bit "0" and the codeword "1" to bit "1". The average number of channel uses per 
bit is thus 1/2 + 1/2-2 = 3/2 . Therefore, the overall number of channel uses necessary for the 
transmission of m bits is upper bounded by + me with arbitrarily small probability for large 
m by the weak law of large numbers (see, e.g., lfT6lO . It follows that an achievable sum-rate is 
given by 

R] + Ro = — i— = — , (5) 
3/2 3 

which is still lower than the upper bound of 1 bit/c.u.. 
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An Optimal Strategy: We now discuss a strategy that achieves the upper bound of 1 bit/c.u.. 
The procedure is based on time-sharing, as driven by the transfer of the energy unit from one 
to the other node. Specifically, each Node j has m bits of information bj t i, 6 j m . If m is large 
enough, from the law of large numbers, approximately half the bits will be zeros and the other 
half will be ones. Since the initial state is Ui ; x = 1, Node 1 is the first to transmit: it sends its 
information bits, starting with bx,i up until the first bit that equals "1". Specifically, assume that 
we have bi t i = = • ■ -fti^-i = and b\^ = 1. Thus, in the iith channel use the energy unit 
is transferred to Node 2. From the (ii + 1)— th channel use, Node 2 then starts sending its first 
bit 6 2 ,i and the following bits until the first bit equal to "1". The process is then repeated. It is 
easy to see that the total time required to finalize this two-way communication is constant and 
equaj^l to 2m and thus the achieved sum-rate is equal to Ri + R 2 = 1 bit/c.u.. 

2) U > 1 Energy Units: In the sum-capacity strategy discussed above with U = 1 energy unit, 
both nodes transmit equiprobable symbols "0" and "1". When there are U > 1 energy units in 
the system, maximizing the sum-capacity generally requires a different approach. Consider the 
scenario with U = 2 energy units: now it can happen that both energy units are available at one 
node, say Node 1. While Node 1 would prefer to transmit equiprobable symbols "0" and "1" in 
order to maximize the information flow to the recipient, one must now also consider the energy 
flow: privileging transmission of a "1" over that of a "0" makes it possible to transfer energy 
to Node 2, leading to a state in which both nodes have energy for the next channel use. This 
might be beneficial in terms of achievable sum-rate. 

Based on this insight, in the following, we propose a coding strategy that employs rate splitting 
and codebook multiplexing. The strategy is a natural extension of the baseline approach discussed 
above for the case U = 1. Each Node j constructs U codebooks, namely Cj\ u , with u G [1,U], 
where codebook Cj\ u is to be used when the Node j has u energy units. Each codebook Cjt u is 
composed of codewords that all have a specific fraction p\\ u of "1" symbols. The main idea is 
that, when the number u of available energy units is large, one might prefer to use a codebook 
with a larger fraction p^ u of "1" symbols in order to facilitate energy transfer. 

2 This equality is approximate if the sequences bj,i, 6 3 , m do not contain exactly half zeros and half ones. 
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Proposition 1. The rate pair (Ri,R 2 ) satisfying 

U 

u=l 

u 

and R 2 < ^2n u H(X 2 \ u ), (6) 

u=l 

where Xj\ u ~ Bern(pjJ, j = 1, 2, /or some probabilities < Piu,P2|« < 1, « = 1 . . . U, wzY/i 
Pi\o = P2\U = 0' i s included in the capacity region C(u). The probabilities n u > 0, it = ... U, 
m © satisfy the fixed-point equations 

K u = 7T«(0O,O|« + 01,l|u) + 7T«-10O,1|« + 7r u+ i0 liQ | u (7) 

wtY/i 7r_i = 7ry +1 = 0, X/«=i = 1» an <^ we ^ ave defined 

0O,O|« = (1 — Pl|u)(l — P2|U-u) 
00,% = (1 -Pl|u)P2|U-u 

0i,o|« = Pi|«(l-Pa[U-u) 
and = Pi|«p 2 [U-u- ( 8 ) 

This proposition is proved by resorting to random coding arguments, whereby codebook Cj\ u is 
generated with independent and identically distributed (i.i.d.) entries Xj\ u distributed as Bern(pj| u ), 
j = 1,2. As introduced above, the idea is that, when the state is = u, Node j transmits 
a symbol from the codebook associated with that state, namely codebook C\\ u for Node 1 and 
codebook C 2 \\j- U f° r Node 2 (which has U — u energy units). Both nodes know the current 
state Ui t i and thus can demultiplex the codebooks at the receiver side. According to the random 
coding argument, the state evolves according to a Markov chain: the system stays in the 
same state u with probability 0o,o|u + 0i,i> (both nodes transmit "0" or "1"), changes to the 
state it + 1 with probability 01^ (Node 1 transmits a "1" and Node 2 a "0") or changes to the 
state u — 1 with probability 0o,i|« (Node 1 transmits a "0" and Node 2 a "1"). The definition of 
the conditional probabilities © reflects the fact that the codebooks are generated independently 
by the two nodes. A full proof is given in Appendix lAl 



October 1, 2012 



DRAFT 



10 



B. Outer Bounds 

In this section, we derive an outer bound to the capacity region C(u). Similar to the standard 
cut-set bound |[T6l Ch. 17], the outer bound differs from the inner bound of Proposition Q] in that 
it allows for a joint distribution <j) Xl X2 \ u of the variables Xu u and X 2 \ u . 

Proposition 2. If the rate pair (Ri, R 2 ) is included in the capacity region C(u), then there exist 

probabilities ir u > with J2 u =i ^ = l > and <t>xiM* - witn Sn,x 2 6{o,i}^i^|« = 1 f or al1 
u G {0, 1, U}, such that 4>i >X2 \o = for x 2 G {0, 1}, X1 uu = for x\ G {0, 1}, condition 

is satisfied, and the following inequalities hold 

u 

R l <Y,^uH{X l \ u \X 2 \ u ) (9) 

u=0 

u 

i? 2 < 5> M # (X 2 | U I X 1K ) (10) 

u=0 

u 

and Rx + R 2 <J2^ H (Xi\ u ,X 2lu ) , (11) 

where variables Xi\ u and X 2 \ u are jointly distributed with distribution 4>xi,x2\u- 

The outer bound above can be interpreted as follows. Suppose that, when the state is = u, 
the nodes were allowed to choose their transmitted symbols according to a joint distribution 
(pxi,x 2 \u — Pr[-X"i,Ai = %i, X 2 ^ k = x 2 ). Note that this is unlike the achievable strategy described in 
the previous section in which the codebook were generated independently. Intuitively, allowing 
for correlated codebooks, leads to a larger achievable rate region, as formalized by Proposition 
2, whose proof can be found in Appendix B. 

C. Numerical Results 

Fig. |2l compares the achievable sum-rate obtained from Proposition 1 and the upper bound 
(fTTI) on the sum-rate obtained from Proposition 2 versus the total number of energy units U. 
As for the achievable sum-rate, we consider both a conventional codebook design in which the 
same probability Pj\ u = 0.5 is used irrespective of the state U\.i = u, and one in which the 
probabilities pj\ u are optimized. It can be seen that using conventional codebooks, which only 
aim at maximizing information flow on a single link, leads to substantial performance loss. 
Instead, the proposed strategy with optimized probabilities pj\ u , which account also for the need 

October 1, 2012 DRAFT 



11 



1.95 
1.9 
1.85 
1.8 

E 

CrT 1-75 
1.7 
1.65 
1.6 
1.55* 
1.50 

23456789 10 
Total energy units 

Figure 2. Achievable sum-rate obtained from Proposition 1 and upper bound i ll It versus the total number of energy units U. 

to manage the energy flow in the two-way communication system, performs close to the upper 
bound. The latter is indeed achieved when U is large enough. 

A remark on the optimal probabilities pj\ u is in order. Due to symmetry, it can be seen that we 
have pi\ u = P2IU-W Moreover, numerical results show that pi\ u increases monotonically as u goes 
from to U, such that p 1 rj > 0.5. In particular, when the number of states U+l is odd, it holds that 
Pi,U/2 — P2.U/2 = 0.5. It is finally noted that the energy neutral transitions (both nodes emitting 
"0" or both emitting "1") occur with equal probability (i.e., (1 — j?i )U )(l — P2,u) = Pi,uP2,u)- 

IV. System Model with Stochastic Replenishments and Losses 

In this section we extend the two-way communication system with energy exchange studied 
above to include energy losses and replenishments, which may occur in different parts of the 
system, as illustrated in Fig. |3] Specifically, the energy units can be lost either while in transit 
through a lossy channel or locally at either node during processing. Similarly, energy units can 
be replenished either by harvesting energy from the channel, e.g., from an interfering signal or 
a source of RF energy, or through a source of power locally connected to the node, e.g., a solar 
panel. All loss and replenishment events are assumed to be independent. As above, we assume 
that the two parties involved have a common clock, and that, at each time, a node can either 
send a "1", which requires one unit of energy cost, or a "0", which does not require any energy 
expenditure. We also assume that Node 1 and Node 2 have energy buffers of capacities B\ and 
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B 2 energy units, respectively, to store the available energy. 



Node 1 Node 2 
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Figure 3. Two-way noisy binary communication with energy exchange. The probabilities of replenishments through the channel 
or locally at the nodes are referred to as pf) , with i ^ j or prs , with i = j, respectively, and similarly for the probabilities 
of losses pfj . See Fig. [4] for an illustration of the channel and Fig. [5] for an illustration of the harvesting process. 




Figure 4. Channel from Node i to Node j. 



Unlike in the previous sections, we assume that the binary channel from Node i to Node j 
with i 7^ j, is noisy as shown in Fig. SI with the probability of of flipping a "0" symbol to 
a "1" symbol and the probability -P^ 1 of flipping symbol "1" to symbol "0". These probabilities 
can be interpreted in terms of replenishments and losses across the channel. To elaborate, let us 
define as the probability of replenishment via harvesting from the channel (for i ^ j), e.g., 
thanks to an RF source that operates on the same bandwidth as the ij link. Moreover, define as 
pfj the probability that an energy unit is lost while in transmit through the channel for the ij 
link. With these definitions, assuming that losses and replenishments are independent, we can 
write the transition probabilities as in Fig. 0] 

Losses and replenishments can also take place locally at the nodes with the probability 
of flipping a "0" symbol to a "1" symbol and the probability P^ of flipping symbol "1" 

(r) 

to symbol "0" at Node i upon reception. Specifically, let us define as p u the probability of 
replenishment at Node i, whereby an energy unit is received by Node i from an external source 
of energy directly connected to the node, such as a solar panel. Note that this energy unit is not 
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received through the channel but is directly stored in the buffer and therefore does not affect the 
decoder, unlike replenishment events over the channel. Moreover, define as pf) the probability 
that an energy unit, while correctly received by the decoder at Node i, is lost during processing 
before reaching the energy buffer. Note that in this case the decoder at Node i correctly records 
a "1", but this energy unit cannot be reused for future channel uses. This event is thus different 
from a loss over the channel in which the decoder at Node i observes a "0" symbol. With these 
definitions, assuming that losses and replenishments are independent and that no more than 
one energy unit can be harvested in each time instant, we can write the transition probabilities 
between the received signal Yj and the harvested energy Hi at Node i as in Fig. \5\ 

— ; o 

N N p(01, =/) W * 

Y H 

1 53 55 " 1 

Figure 5. Statistical relationship between the received signal Yi and the energy Hi harvested by Node i. 

Based on the discussion above, at any given time instant k, with k G [l,n], the state of the 
system {U\^i U 2 ,k) G N 2 is given by the current energy levels U\ t k and U 2 ,k in the buffers of Node 
1 and Node 2, respectively. By the capacity limitations of the buffers, we have the inequalities 
Mi G [0,Bi] and u 2 G [0,B 2 ] for each channel use k G [l,n]. The transmitted symbols are 
limited as per (OQ). 

The channel is noisy with transition probabilities as in Fig. [3] Moreover, the relationship 
between received signal and harvested energy is as in Fig. |5] Therefore, the state of battery at 
Node 1 for k G [1, n] evolves as follows 

Ui t k = (U^k-i — ^i,fc-i) + H ik _i. (12) 

Similar to Sec. II, we use the following definitions for an (n, R 1 , R 2 , B 1 , B 2 ). Specifically, 
the code is defined by: the buffer capacities Bi and B 2 ; two sequences of encoding functions, 
fi t k for k G [1, n] and i = 1,2, which map the message Mi, the past received symbols Y^~ x 
along with the past and current states (U±, C/|) into the currently transmitted signal X itk G Xu ik \ 
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two decoding functions for % = 1,2, which map the received signal Y™, the local message 
Mi and the sequence of states 11^,1/ 2 into an estimate Mj of message Mj of the other node 
j ^ i. Achievability is defined as in Sec. II. Finally, the closure of the set of all the rate pairs 
(Ri, R 2 ) is defined as the capacity region C (-Bi,£>i), where subscript "0" denotes the fact that 
the capacity region generally depends on the initial state. 

Remark 3. In the definition of code given above, we have assumed that the nodes can track 
the state of the batteries U 2 ,k) at both nodes. We refer to this scenario as having Global 

Energy Information (GEI). We remark that in the presence of losses and replenishment, the 
nodes generally cannot track the amount of energy available at the other node based only on the 
knowledge of the received signal. Instead, information about the state of the other node needs 
to be acquired through additional resources such as control channels or appropriate sensors. In 
general, the assumed model with GEI can thus be thought of providing a best-case bound on 
system performance. In Sec. IV-Cl we will study the scenario, referred to as having Local Energy 
Information (LEI), in which each node is only aware of the energy available in its own local 
battery. 

V. Generalizing the Inner and Outer Bounds 

In this section, we first propose a communication strategy that leads to an achievable rate region 
by generalizing the approach discussed in Sec. HID The outer bound of Sec. HH] is similarly 
extended. While the strategy at hand is based on GEI (see Remark [3]), we then discuss an 
achievable strategy with LEI in Sec. IV-Cl and present some numerical results in Sec. IV-Dl 

A. Transmission Strategy 

The proposed strategy is an extension of the approach put forth in Sec. HU] and operates as 
follows. Node i, with i — 1,2, constructs one independent codebook CiU ul , U2 ) f° r eacn possible 
sate (ui, U2) £ [0, B{\ x [0, B 2 }. As in Sec. Hill at each time k, if the state is (L^, U 2 ,k) = (wi, u 2 ), 
then Node i transmits the next symbol from the codebook CiU ul>U2 )- At the end of the last channel 
use, each node, being aware of the sequences of states, can demultiplex the transmission of the 
other node and decode the messages encoded in all the (B 1 + 1)(B 2 + 1) codebooks. 

The codebook of Node i corresponding to state (ui,u 2 ) is generated by drawing each bit 
independently with a given probability Pi\( Ul ,u 2 ) f° r i = 1, 2 and all states (ui,u 2 ) £ [0, Bi] x 
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[0, B 2 \. Note that, due to CQ), we have Pi|(o,« 2 ) = for all u 2 G [0,B 2 ] since, when Ui tk = 0, 
Node 1 has no energy available and thus must transmit a "0" symbol; and similarly we have 
P2|(iti,o) = f° r ai l u i e [0,-Bi]. Given the probabilities Piu ui ,u 2 ) f° r i = 1, 2 and all states 
w 2 ) G [0, B x ] x [0, £ 2 ], the (B 1 + 1)(B 2 + 1) x (S X + 1)(S 2 + 1) transition probability matrix P 
can be obtained that contains the transition probabilities from any state (u±, w 2 ) G [0, B]\ x [0, £? 2 ] 
to any state (u'^u^) G [0,£>i] x [0,5 2 ]. These transition probabilities depend on the parameters 
(Pij iPij)> (Pii\Pii)> and PiKu u u 2 ) for h3 = !)2 and (ui,w 2 ) G [0, Si] x [0,S 2 ], as detailed in 
Appendix O 

5. Inner and Outer Bounds 

In order to derive the rates achievable with this strategy, denote as 7T( U1>U2 ) the average fraction 
of channel uses k such that we have (Ui^, ^2,fe) — (^1,^2) for all states («x,u 2 ) G [0,i?i] x 
[0,S 2 ], as done in Sec. [Till Note that J2( Ul « 2 ) 7r ( M i. M 2) = 1- This function is also referred to as 
the steady-state probability and can be calculated as the limit 



where 7T( UljU2 ) is the {B\ + 1)(-B 2 + 1) x 1 vector containing the steady-state probabilities 7T( UliU2 ) 
for all states (ui,u 2 ) G [0, Bi] x [0, B 2 ) and we recall that P is the transition probability matrix. 
Vector 7r(l) accounts for the initial state and is thus a vector of all zeros except for a one in 
the entry corresponding to the initial state. We note that the limit in (TT3T) always exists for the 
model studied in Sec. [HI] (for all non-trivial transmission probabilities), and is given by CZ])-®. 
The same is generally true here apart from degenerate cases. However, the transition matrix (TT3T) 
is possibly reducible, and thus the calculation of the limit generally requires the factorization of 
the matrix according to the canonical form for reducible matrices. We refer to lfT8l ch. 8] for a 
detailed discussion on the existence and calculation of the limit (TT3T) . 

Proposition 4. Assuming that the limit rfiil) exists, the rate pair {R\, R 2 ) satisfying the inequal- 
ities 



(wii«2) 



lim P fc 7r(l) 



(13) 





(«i,«2)e 

[0,B 1 ]x[0,S 2 ] 



and R 2 < 





(14) 



(«i,«2)e 

[0,B 1 ]x[0,S 2 ] 
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for some transmission probabilities Pi\( ui ,u 2 )> f or i = 1,2 and (ui,u 2 ) G [0,i?i] x [0,_B 2 ] is 
achievable, where we have denoted as Xi\i UljU2 -\ as the Bernoulli variable Bern(pi\i Uim \). We 
also have 

I(Xi\(u lt u 2 )] Y 2 ) = H ((1 - pi|( Ul>U2 ))P i + Pi|(ui,« 2 )-Pii) 

- [p 1](ui U2) H{P 11 ) + (1 - p ll{ui U2) )H{P 01 )] (15) 

and similarly for I(X 2 ^ UliU2 y, Y"i). 

Remark 5. The achievability of the rates in ([14)) can be proved by adopting the multiplexing 
strategy described above and following the same main steps as in Appendix A. Here, we also 
point out that the achievability of (fl4"l) under the assumption that the limit (fl"3T) exists is a direct 
consequence of [f9l Lemma 12.3.1]. 

An outer bound can be also derived by generalizing Proposition 2. In particular, following 
similar steps as in Appendix B, one can prove that an outer bound is obtained by allowing for 
joint probabilities, rather than product distributions as in Proposition 1. Moreover, one can add 
the sum-rate constraint that generalizes (fTTj) as 

Ri + R 2 < ^2 n (ui,u2)I(.Xi\(u u u 2 ), X 2 \( UUU2 y, Y 2 ), (16) 
{ui,u 2 )e 

[0,B 1 ]X[0,B 2 ] 

where Xi|( Uli „ 2 ), X 2 \( Ul ,u 2 ) are jointly distributed. 

C. Local Energy Information 

In the discussion above, we have assumed GEI, that is, each node knows the full current 
energy state (Ui^, U 2) k) (see Remark [3]). In this section, we consider instead the scenario with 
LEI, in which Node 1 only knows its local energy level U\ and Node 2 only knows U 2 . 

We first observe that the energy can be considered to be the state of the link 12 at channel 
use k, since it affects the available input symbols via © (and similarly for U 2 ^ and link 21). 
Therefore, the model at hand falls in the category of channels with states in which the state is 
known only at the transmitter. For these channels, under the assumption that the state sequence is 
i.i.d. and independent of the transmitted signal, it is known that so called Shannon strategies are 
optimal [16, Ch. 7]. In the model under study, unlike the conventional setting, the state sequence 
U'i (and U 2 ) is neither i.i.d. nor independent of the transmitted signal X" (and X%). Therefore, 
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Shannon strategies are generally not optimal. We will see below that they can be nevertheless 
used to lead to non-trivial achievable rates. 

Following Shannon strategies, we draw auxiliary codebooks made of independent and i.i.d. 
codewords V™ and V 2 n using pmfs p{y\) and p(v 2 ), respectively. Each symbol V, fc for Node j 
and time instant A; is a vector consisting of Bj bits. The main idea is that, at each time k, Node 
j transmits the bit in Vj f. corresponding to the current state Uj^. Note that the latter can take 
Bj possible values at which the transmitted signal is non-trivial (for Uj^ = 0, we necessarily 
have Xj jk = 0). 

At the receiver side, the decoder at Node 2 uses joint typicality decoding with respect to the 
distribution p{v 1,1/2), which is given as 



where tt(ui) is the marginal distribution of the steady-state probability of the Markov chain 
induced by the random coding strategy and the evolution of the system, as discussed above (see 
also Appendix C). Following standard information theoretic considerations, we obtain that the 
rate pair (Ri,R 2 ) satisfying 



for some pmfs p(v 1), p(v 2 ) is achievable, where p(v 1, y 2 ) is as in (fT7l) and similarly for p(v 2 , yi). 
Regarding the details of the proof, being based on conventional tools (see [16, Ch. 3]), here we 
simply point out that it is based on the ergodicity of the Markov chain, which allows to conclude 
that the error event in which the correct codeword is not jointly typical takes place with negligible 
probability; and the packing lemma in IfTBl Lemma 3.1], which entails that the error events due 
to mistaking other codewords for the correct one have also negligible probability^. 

Remark 6. In the strategy proposed above, each node adapts the choice of the current transmitted 
symbol only to the current local energy state. A potentially better approach would be to perform 
adaptation based on a local state that includes also a number of past energy states of the node, 

3 The packing lemma does not assume that the received signal be i.i.d. and thus applies to our scenario (see 1161 Lemma 3.1]). 




(17) 



(18a) 



andi? 2 < /(V^n) 



(18b) 
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along with the current one, and/or current and past received signals. This aspect is not further 
explored in this paper. 

D. Numerical Results 

In this section, we present some numerical examples in order to assess the impact of replen- 
ishment and loss processes. We assume that each node has the ability to store only one unit 

(r) (r) 

of energy i.e., B\ = B 2 = 1, and we consider a symmetric system with p 12 = p 2 { = Pr,c, 
Pii = P22 = Pr,n, P12 — P21 = Pi,c an d P11 = P21 — Pl,n> where the subscripts "c" and "n" stand 
for "channel" and "node", so that, e.g., p r n is the probability of replenishment locally at a node. 
We first assume GEI. 

Fig. 5 shows the sum-rate obtained by summing the right-hand sides of (fl4l) . optimized over 
the probabilities pi\( Ul u 2 ) and p2\( Ul u 2 ) for ai l states (u\, u 2 ) G [0, B\] x [0, B 2 ] versus the replen- 
ishment probability on the channel p rc (see Fig.H]) for two cases, namely p rn = 0,Pi, n ,Pi,c — 0.1 
and p r>n = 0,Pi, n ,Pi,c = 0.3. We also show in the same figure the steady-state probability n^i) 
of state (ui,u 2 ) = (1, 1) corresponding to the optimal values of Pi\( Ul u 2 ) and p 2 \( Ul u 2 )- It is seen 
that increasing the probability p r c increases the chance of being in state (ui,u 2 ) = (1, 1), due 
to the increased availability of energy. However, increasing p r c has also the deleterious effect 
of flipping bits on the channel from "0"s to "l"s with larger probability. It is seen that, in the 
regime in which p rc is sufficiently small, and the system is energy-limited, increasing p r c is 
beneficial, while for p r c large enough the second effect dominates and the achievable sum-rate 
decreases. 

We now turn to assessing the effect of local replenishment at the node. Specifically, Fig. 6 
shows the optimized sum-rate versus p r , n . As we can see, increasing p rri improves the sum-rate, 
since it enhances the probability of being in state (ux, u 2 ) = (1, 1), without any side effect since 
it does not impair the channels. 

Fig.[8]and Fig.[9]show the effect of loss events on the channel and at the nodes, respectively. We 
show both the sum-rate and the optimal transmission probability pi\(\ : o), which equals the optimal 
probability p 2 |(o,i) by symmetry. The latter is also compared with the transmission probability that 
maximizes the mutual information I(Xiu 0>1 y, Y 2 ) in ([141) and that is thus capacity achieving. It is 
noted that this is the probability that maximizes the information rate when there are no energy 
limitations. As it can be seen, by comparing Fig. [8] and Fig. HI increasing the loss probability 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

Pr.C 



Figure 6. Sum-rate R sum and steady-state probability nv^i) versus the probability p nc of replenishment on the channel (see 
Fig. H. 





0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

Pr,n 



Figure 7. Sum-rate R SU m and steady-state probability ttm i) versus the probability p r , n of replenishment at the node (see Fig. 

H). 

both on the channel and at the node decreases the sum-rate, although the rate of this decrease 
is larger for the latter, since, similar to the discussion above, a loss at the node does not affect 
the channel. Moreover, for small pi >c and pi <n , the transmission probability pi|(i,o) is close to the 
capacity-achieving probability, while for larger loss probabilities p lc and pi >n , it becomes smaller 
than the capacity-achieving probability. 

We now consider the effect of LEI. Fig. [TO] compares the sum-rate achieved with GEI and 
LEI versus the replenishment probability p TjC on the channel. As it can be seen, LEI entails 
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Pl,c 



Figure 8. Sum-rate R S um and optimum transmission probability Pik^o) versus the probability p;, c of loss on the channel. 
Dotted lines show the capacity achieving probability. 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



2 0.45 

or 0.4 









Cap. ach. p lc = 0,p, c =0.1 

Cap.ach. p u = 0,p r , = 0.3 







0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Figure 9. Sum-rate R S um and optimum transmission probability Pi|(i,o) versus the probability pi, n of loss at the node. Dotted 
lines show the capacity achieving probability. 



a significant performance loss with respect to GEL To gain some insight as to the reasons of 
this loss, the figure also shows the optimal transmission probabilities Pi|(i,o)> Pi|(i,i) w i m GEI 
and the probability p(v\) = pi\\, that is the probability of transmitting "1" if the local battery 
contains energy, for LEI (T4 is a Bernoulli variable since B x = B 2 = 1). With GEI, the nodes 
can adapt the transmission strategy to the energy state of both nodes and thus choose different 
probabilities Pi|(i,o) and Pi|(i,i), while with LEI the nodes are forced to choose a single probability 
irrespective of the state of the battery at the other node. 
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0.35 




Figure 10. Sum-rate R aU m and optimum transmission probabilities for Global Energy Information (GEI) and Local Energy 
Information (LEI) for p ijC = 0.3 and p r , n = p i>n = 0. 



VI. Conclusions 

Energy and information content are two contrasting criteria in the design of a communication 
signal. In a number of emerging and envisaged communication networks, the participating nodes 
are able to reuse part of the energy in the received signal for future communication tasks. 
Therefore, it becomes critical to develop models and theoretical insights into the involved trade- 
offs between energy and information exchange at a system level. In this work, we have taken 
a first step in this direction, by considering a two-way channel under a simple binary "on-off" 
signaling model. The derived inner and outer bounds shed light into promising transmission 
strategies that adapt to the current energy state. It is emphasized that conventional strategies 
based solely on the maximization of the information flow entail substantial losses. 

The results presented in this paper call for further studies on different fronts. One is the 
development of better models which strike a good balance between adherence to reality and 
analytical tractability. A second is the development of better communication strategies for the 
practical scenario in which the energy state of the network is not fully known at the nodes. 

Appendix A 
Proof of Proposition 1 

1) Code construction: We generate U codebooks for each Node j = 1,2, namely Cj\ u , with 
u G [1,U]. The codebook Cj\ u for u > has Kj )U codewords, each consisting of n^ u symbols 
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%j,u,i £ {0)1}' which are randomly and independently generated as Bern(pj| u ) variables, with 
I — 1, 2, n 3 - jU and n JiU = for some < Sj :U < 1. We denote the codewords as x^ u (mj jU ) 
with m JiU G [1, Kj iU ]. Note that the parameter Sj iU does not depend on n, and hence, if n — > oo, 
then we have n jjU — >■ oo for all j, it. We set 2 nR j = n«=i -^?>' while the relations among the 
remaining parameters (Kj jU ,5j :U ,pj\ u ) will be specified below. 

2) Encoding: Each node performs rate splitting. Namely, given a message Mj G [1,2"^'], 
Node j finds a U— tuple [ra^i, ra^u] with m J)U G [1, that uniquely represents Mj. This is 
always possible since we have 2 nR j = Yl u =i Kj,u- Then, the selected codewords x^ u u (rrij tU ) for 
it G [l,u] are transmitted via multiplexing based on the current available energy. Specifically, 
each Node j initializes U pointers Z^i = lj >2 = ■ ■ ■ = Z ■ y = 1 that keep track of the number of 
symbols already sent from codewords x^^irrij^), x 1 ^ '(rrijp), x. 3 ^(m^\j), respectively. At 
channel use i, if the state is U\ :i = u, then the nodes operate as follows. 

• Node 1: If u — 0, then xi^ = 0. Else, if h iU < ni jU , Node 1 transmits x\ yi = xi jUi i lu (mi jU ) 
and increments the pointer l ljU by 1. Finally, if l ljU = n 1:U + 1 the pointer v\ jU is not 
incremented, and the transmitter uses random padding, i.e., it sends x lti = 1 with probability 
Pi yU and xi :i = otherwise. 

• Node 2: If u — U (i.e., no energy is available at Node 2), then x 2ji = 0. Else, if Z 2 ,U-u < 
"2,U-«. Node 2 transmits x 2 ,i = x 2 u 2l (m 2 n_J and increments the pointer Z 2 n- U 

2, U — it 

by 1. Finally, if Z 2 ,U-« = n 2,U-u + 1, the pointer Z 2 ,U-u i s n °t incremented, and Node 2 

sends x 2ji = 1 with probability p 2j u-« and x 2 ,i = otherwise. 
The random padding method used above is done for technical reasons that will be clarified 
below. 

3) Decoding: We first describe the decoding strategy for Node 2. By construction, the nodes 
are aware of the state sequence and thus can determine the ordered set 

K = {i\Ui,i = u}, (19) 

of channel use indices in which the state is u with u G [0,U]. For all u G [1,U], if \M U \ > n 1:U , 
then Node 2 takes the first n\ :U indices i Ujl < i Uj2 < • • • < i u ,n 1>u from the set Af u and obtains 
the list of messages toi iU G [1,Ki )U ] that satisfy 5i, u ,fe(mi )U ) = x^ iuk for all k G [l,ni )U ]. Note 
that the list cannot be empty due to the fact that the channel is noiseless. However, it contains 
more than one message, or if \M a \ < n 1:U , then Node 2 puts out the estimate rhi :U = 1. Instead, 
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if the list contains only one message m 1(1 , then Node 2 sets rh 1>u = m lu . Finally, the message 
estimate is obtained as rh\ = [m^, ...,rh x xj\. 

Node 1 operates in the same way, with the only caveat that the nth codebook C 2 \ u of Node 2 
is observed at channel uses in the set Njj- u for u E [1, U]. 

4) Analysis: We evaluate the probability of error on average over the messages and the 
generation of the codebooks, following the random coding principle. From the definition of 
the decoders given above, the event that any of the decoders is in error is included in the set 
£ = Ui=i^US=i(4« U Efl), where: (f) Sf} is the event that \M U \ < n ljU for j = 1 and 
that |jVg_J < n 2jU for j = 2, that is, that the number of channel uses in which the system 
resides in the state in which the codeword x^ u u (rrij !U ) from the codebook Cj tU is sent is not 
sufficient to transmit the codeword in full; (if) Efl is the event that two different messages 
m' ju} m" u E [l,K j:U ] are represented by the same codewords, i.e., x^ u (m' ju ) = ^"(m^). 

The probability of error can thus be upper bounded as 

2 U 
j=l u=l 

In the following, we evaluate upper bounds on this terms. 

It immediately follows from the packing lemma of [16] that Pr[£-^] — > as rij, a — > oo as 
long as 

l0g2 ^ u <H( Pj]u )-5(e) (21) 



with 5(e) — > as e — > 0. For analysis of the probabilities Pr[£j^], we observe that, under 
the probability measure induced by the described random codes, the evolution of the state Un 
across the channel uses i E [1, n\ is a Markov chain with U + 1 states. Specifically, the chain is a 
birth-death process, since, if the state is XJ\^ = u in channel use i, the next state U\^ + \ can only 
be either u — 1 or u+ 1. More precisely, let q u \ w = Pr(L r i ii+ i = u\Ui yi = w) be the transition 
probability. Note that, due to the use of random padding, the transition probability q u \ w remains 
constant during all n channel uses, so that the Markov chain is time-invariant. 

We now elaborate on the Markov chain Ui ti . To this end, we first define as (j) Xl ,x 2 \u, where 
Xi,x 2 E {0, 1} be the joint probability that Node 1 transmits X 1A = x% and Node 2 transmits 
X 2 ,i = x 2 during the ith channel use in which the state is Ui ti = u. We can now write the 
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non-zero values of the transition probability q u \ w as follows: 

Qu,u-1 — 01,O|u Qu,u+1 = 0O,l|u 

Qu,u = 1 — Qu,u-l ~~ Qu,u+i (22) 

With a slight abuse of the notation and noting that 01,010 = 0i,i|o = an d 0o,i|U = 0i,i|U = 
the expressions above also represent the transitions for the two extremal states u = and u = U, 
as they imply g |-i = and gij|U+i = 0- 

If Pi,o = P2,o = and < pi >u ,P2,u < 1 for all u > 0, then it can be seen that the Markov 
chain is aperiodic and irreducible, and thus there exist a unique set of stationary probabilities 
7To, iti, • ■ ■ , 7TQ, which are given by solving the linear system, defined by taking U equations of 
type © for u = . . . U — 1 and adding the condition Y^=o = 1- 

We are now interested in the statistical properties of the set \M U \ of channel uses in which 
the state satisfies U\ = u. Using the ergodic theorem and the strong law of large numbers ifTTl 
Theorem 1.10.2], it can be shown that limn^oo = tt u with probability 1. Therefore, if we 
choose: 

h, u = k,\J-u = n{ir u - e) (23) 

then Pr[£^2] = Prff^u^J can be made arbitrarily close to as n — > oo. This concludes the 
proof. 
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Appendix B 
Proof of Proposition 2 

Consider any (n, Ri, R 2 , U) code with zero probability of error, as per our definition of 
achievability in Sec. HB We have the following inequalities: 

nR x = H(M 1 ) = H{M X \M 2 , U 1A = u 1A ) 
( = } H(M 1 ,X?,U?\M 2 ,U l!l = u 1 , l ) 
^H(X?,U?\M 2 ,U ltl = Ul!l ) 

n 

= H(X lti , Ut^Xt 1 , Ut\ M 2 , U ltl = ui tl ) 

n 

= £ H{U ltt \X{-\ U{-\ M 2 , U 1A = u lt i) 

i=l 

+ H(X u \X{- 1 ,Ui,M 2 ,U lil = u ljl ) 

n 

( = } H(X 1 , i \Xt 1 , U[, M 2 , U ltl = «x,i) 
i=i 

< ^^(X^C/^X^) 

i=l 

{ ^ H{X 1 \U l ,X 2 ,Q) 

<H(X 1 \U 1 ,X 2 ), (24) 

where (a) follows since X™, [7™ are functions of Mi, M 2 and u^i, (b) follows since H(M± (Xj 1 , [7™, 
M 2 , U11 = = holds due to the constraint of zero probability of error; (c) follows since 
Ui t i is a function of X{ _1 , M 2 and ui t i*, (d) follows by conditioning reduces entropy; (e) follows 
by defining a variable Q uniformly distributed in the set [l,n] and independent of all other 
variables, along with X x = Xiq, X 2 = X 2 q and Ui = Uiq. 

Similar for nR 2 we obtain the bound nRi < H(Xi\Ui,X 2 ). Moreover, for the sum-rate, 
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similar steps lead to 

n{R x + R 2 ) = H(Mi, M 2 ) = H(M 1 , M 2 |C/ M = u 1A ) 
= H{M l M 2 ,X^X^U r ^\U ltl = u hl ) 
= H{X?,X$,U?\U ltl = u ltl ) 

n 

= H(U lti \Xi-\ Xt 1 , U{-\ M 2 , U 1A = Ul)1 ) 
i=i 

+ H(X 1>h X 2 ,\X\-\X l 2 -\ Ul M 2 , U ltl = u ltl ) 

n 

>^ff(X M ,X 2ii |C/ M ) 

i=l 

= H(X U X 2 \U 1 ). (25) 

Let us now define n u = Pr[C/i = u] and (f> xl>X2 \ n = Pr[Xi = x±,X 2 = x 2 \Ui = u] for 
i, j E {0, 1} and for all U\ E {0, 1, U}. Probability conservation implies that the relationship 
CO) be satisfied. This concludes the proof. 

Appendix C 

Transition Probabilities for the Model in Sec.[V] 
Here we discuss the transition probability matrix P used in Sec. |V} To this end, define as 
Qij b ^ (hJ — 1) 2) the probability that Hj = a E {0, 1} energy units are added to the battery at 
Node j conditioned on Node i sending symbol Xj = b E {0,1}, for i ^ j, namely 

(26a) 
(26b) 
(26c) 
(26d) 

Note that these transition probabilities correspond to the cascade of the channels in Fig @] and 
Fig [51 Based on these probabilities, we can now evaluate all the possible transition probabilities 
from state (ui, u 2 ) to any other state (u' x , u' 2 ). We start with u\ E [1, B 1 — 1] and u 2 E [1, B 2 — 1] 
for Bi, B 2 > 1 whose outgoing transition probabilities are illustrated in Fig.QT] The "boundary" 
states with Uj = or Uj = Bj for some j = 1, 2 are discussed later. 
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Figure 11. The outgoing transition probabilities from a state (u\,u 2 ). 

By the stated assumptions, the state (ui,u 2 ) can only transit to state (ui + ii,u 2 + 12) with 
^1^2 £ { — 1?0, 1}, so that the energy in the battery is increased or decreased by at most one 
energy unit. Therefore, for the "non-boundary" states (ui,u 2 ) with u\ G — 1] and u 2 E 

[1, B 2 — 1], the probabilities in Fig. [TT|can be easily obtained as 

Po,o{ui,U 2 ) = Pl\{u 1 ,u 2 )P2\{u 1 ,u 2 ){Q 2 l ) Ql 2 l) ) + Pl\{u 1 ,u 2 )p2\{u 1 ,u 2 ){Q 2 l ) Ql2 ] ) 

+Pl\(u 1 ,u 2 )P2\(u 1 ,u 2 )[Q 2 l Q\ 2 ) + Pl\(u 1 ,u 2 )P2\(u 1 ,u 2 ){Q21 Ql2 J) ( 2/ ) 

P +t0 {Ul,U 2 ) = Pl|( UllU2 )P2|(t 11 ,t 12 )(Q21 1)( 5l 2 1) ) +Pl\(u 1 ,u 2 )P2\(u 1 ,u 2 ){Q 2 l ) Q i l 2 ) ), (28) 

Pq,+\UI,U 2 ) = Pl|(«i,«2)P2|(« 1 ,« 2 )lV21 V12 J + Pl\(u u u 2 )P2\(u u u 2 )W21 Vl2 J) ( 2y ) 

■P+,+ («l,«2) = Pl|(ui,« 2 )^|(ui,«a)(Qsil V12 ), (30) 

P+_(mi,M 2 ) = Pl|(«i,«2)P2[(«i, U2 )( ( 521 1)( 5l2 0) ) ) (31) 

P-,+ («l,«2) = Pl|(«i,«2)P2[(ui lU2 )(Q21 0)( 5l2 )> (32) 
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Po,-(Ul,U 2 ) = Pl|(«i,«2)P2|(«i,«a)(Q21 Ql2 ) +Pl|(«i,«a)i ? 2|(«i,«2)(Q21 Ql2 )> ( 33 ) 

^-,0(«l,«2j = Pl|(«i,«2)P2|(« 1 ,« 2 )W21 V12 J +Pl|(ui,«2)P2|(ui,u 2 )lV21 Vl2 )■> ( 34 ) 

and 

P__(m1,U 2 ) = Pl\(u l ,u 2 )P2\(ui,u 2 )(Q21^Ql2 )j (35) 

where we recall that Pi|( Ul ,u 2 ) is the probability of sending a "1" symbol by Node z given the 
state («i,« 2 ) and Pi|( Wl)Wa ) = 1 -Pi|( UljU2 )- 

For the "boundary" states (/Ui,^) with «i and/or u 2 equal to the outgoing transitions in 
Fig. QT] and probabilities in (l271)-d35l) still hold since the transitions to states with energy less 
than zero are dis enabled by the conditions Pi\( Ul ,u 2 ) — for Uj = 0, j = 1,2. Instead, if 
Mi = B\ and m 2 G [0, 5 2 — 1] then the probabilities in (l33l)-(T35l) remain the same, but we 
have P +; _ = P +)0 = P+.+ = and P .o equals the sum of the right-hand sides of (127|) and 
(T28T) . while P + equals the sum of the right-hand sides of (T29| ) and (l30l) . and P _ equals the 
sum of the right-hand sides of (|33l) and (I3TI) . The transition probabilities from the states with 
U\ G [0, Bi — 1] and u 2 = B 2 follow in a symmetric fashion. Finally if u\ = Bi and u 2 = B 2 , 
then we have P_ + = P + _ = P +j0 = P ,+ = P+,+ = and P 0i0 is the sum of ([27]), (EHJ), 
and ([29]). while P _ is the sum of ([33]) and (HB and P_ j0 is the sum of ([34]) and (|32[) . 

By using the transition probabilities defined above, one can easily construct the transition 
matrix P of the corresponding Markov chain. For instance, for the case B = 1, we can write 
the transition matrix as 

P ,o(0,0) P ,-(0,1) P_ )0 (1 ) 0)P_ ) _(1,1) 

P o , + (0, 0) P ,o(0, 1) P_ >+ (1, 0) P_, (l, 1) 
P= , (36) 

P +i o(0,0)P +) _(0,l) P ,o(l,0) P ,-(l,l) 

P +)+ (0,0)P + , (0,l) P ,+(1,0) P 0)0 (l,l). 

where the column index represents the the initial state and the row index the final state. 
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